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1.0 SUMMARY 

An examination of the factors which modify the simulation of a 
constraint in the motion of the aft attach points of the Orbiter 
and External Tank during separation has been made. The factors 
considered were both internal (spring and damper constants) and 
external (friction coefficient and dynamic pressure). The 
results of this study show that an acceptable choice of spring/ 
damper constant combinations exist over the expected rang? of the 
external factors and that the choice is consistent with a prac- 
tical integration interval. This constraint model is also shown 
to produce about a 10% increase in the relative body pitch angles 
over the unconstrained case whereas the MDC-STL constraint model is 
shown to produce about a 38% increase. 


2.0 INTRODUCTION 

A constraint model was developed by TRW for use in the Space Vehicle 
Dynamics Simulation (SVDS) program to provide a simulation of the 
binding characteristics of the ball and socket associated with the 
aft attach points of the Orbi ter/External Tank during separation. 

This model, which is called ABIND, uses a spring-damper-niass 
system to model the force and moment interactions in the ball 
and socket. These interactions, which occur between the separat- 
ing vehicles after release of the restraining mechanisms but be- 
fore actual loss of physical contact must then be incorporated in 
to the equations of motion of each of the two bodies. As in all 
spring-damper-mass systems the force histories transmitted to the body 
masses are a function of the spring stiffness and damping used in the 
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simulation and the choice of spring and damper constants becomes 
representative of the quality of the simulation. This study 
will show the relationship between the choice of spring/ 
damper constants and the corresponding responses to them for a range of 
probable operating conditions. The primary operating conditions 
include the friction coefficient between the aft ball and socket 
and the dynamic pressure at separation. A comparison with an 
unconstrained simulation and an alternate constraint model are 
also briefly examined. 

3.0 DISCUSSION 

The ABIND constraint model provides a simulation of the compon- 
ents of the normal and frictional forces and torques that are 
developed within the ball and socket of the two aft attach 
points. This simulation uses a spring-damper-mass system to 
provide the force and torque components, but a tuning of the 
system is required to select the proper spring and damper 
constants. This selection must consider the force and moment 
profiles generated by the motion, the integration interval 
required by the computer to convert one force cycle 
into motion, and the resultant motion of the center of the 
ball as it leaves the socket. 

Integration Interval - The aft attach of the Orbiter and External 
Tank is modeled by a simple parallel spring-damper-mass system 
and is, therefore, amenable to a theoretical estimate of the 
Integration interval. From the basic formulation of spring-mass 
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systems, the rotational frequency of the oscillating forces can 
be derived as: 


% 


where 


p = f2Kg(W^ - W 2 ) (rad/sec) 

\ W^W2 

= Weight of Arbiter (lb.) 

= Weight of External Tank (lb.) 

K = Spring Constant (Ib/ft.) 

2 

g = Gravitational Acceleration (ft/sec ) 


The period (T) can then be found as: 



(sec. ) 


The minimum number of points required to define the cyclic 
forces in a 4th order Runge-Kutta integration is estimated at 
six per period. Thus, the maximum integration interval (At) 
is then found to be: 

T TT / W,W, 

" 6 " aV 2Kg(W^ - W 2 ) 


As shown in Figure 1 a spring constant no greater than 

2.08 X 10^ Ib/ft. can be used with a maximum integration interval 

of .01 seconds for tank weightj exceeding 80,000 lbs. 

As an example, a safe comparison of an actual spring constant- 
integration interval relationship was explored by comparing the 
accuracy of a At of .01 and .002 seconds for a typical spring 
constant of 10^ Ib/ft. These results, which are presented in 
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Figure 2 shows that there is a very close match in the torques 
of the External Tank during separation. Even better results 
are obtained by comparing the integrated trajectory of the 
initial x-z motion of the ball in the socket, as in Figure 3. 
Thus, an integration interval of .01 seconds is capable of 
accurately following the dynamics corresponding to a spring 
constant of 10^ Ib/ft. 

T ypical Data - Having selected an integration interval of .91 
seconds, a typical data set can now be provided to describe the 
force, torque and trajectory outputs for an arbitrary input. 

This input uses the 10® Ib/ft. spring constant along with a 
damper -constant (C) of 10® Ib-sec/ft, a friction coefficient (y) 
of .2 and a dynamic pressure (q) of 10 Ibs/ft . Other intial 
conditions and mass properties which are used throughout the 
study are contained in Table 1. 

Figure 4 shows the Orbiter reaction force profiles for both 
the left and right attach points. It is apparent that the 
out-of-plane initial conditions cause the left attach to 
separate about .2 seconds before the right. It also shows that 
the axial force is far larger than either the vertical or 
lateral forces and that none of the force histories have any 
significant oscillations. 

The torques about the x, y, z axes thru the Orbiter center-of- 
gravity are shown in Figure 5. This figure shows that the 
forces acting about the Orbiter c.g. produce significant torques 
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ABIND ANALYSIS 

INTIAL CONDITIONS AND MASS PROPERTIES FOR STUDY 
INITIAL CONDITIONS* 


= -4 deg. 

= 6^ = -0.5 deg. 

0^ = -0.25 deg/sec. 

4^. = 0.5 deg/sec. 

= -0.5 deg/sec 


14 psf 
195,927 ft. 


Vg = 6750 ft/sec. 

ipg = 90.0 deg. 

q = 2 psf 10 psf 

h = 245,105 ft. 204,970 ft. 


MASS PROPERTIES 


ORBITER 

TANK 

Weight (lb.) 


230942 

80000 

CG Location (Inches) 

(in Vehicle Reference 

X 

V 

1111.9 
n ^ 

1363.2 

2.0 

Coordinate System) 

Y 

U.o 

z 

380.6 

418.5 

Inertias (10® slug-ft^) 

^xx 

0.8378 

0.396 


*yy 

6.489 

4.588 


6.70 

4.556 


■xy 

■xz 

0.003635 

0.0 


0.1615 

0.0 

m 


0.002196 

0.0 


♦These initial conditions are identical with those provided MDTSCO 
by MDAC-STL for out-of-plane case JJ and IV dated 24 August 1974. 
Constraint model was used only for case IV. 
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about the y and z axes. The torque profile about the y-axis 
(pitching torque) Is slmliar In shape to the axial force profile 
although It also contains components of the vertical force In 
it. The torques about the z-axis (yawing torque) for each of 
the attach points are similar in shape but of opposite sign and 
hence tend to cancel each other whereas the pitching torques are 
additive. The torques about the x-ax1s (rolling torque) are 
very small by comparison. It can thus be shown that the single 
most significant reaction parameter for assessing spring/damper 
effectiveness can be demonstrated by the pitching torques. 

The trajectories generated by the integration of the previously 
described forces and torques profiles are also indicative of 
the spring/damper constant response. Two types of trajectories 
are presented in Figure 6. The first shows the translation of 
the ball from the socket and the second shows the rotation of 
the ball as it is displaced vertically. These plots show that 
the left and right ball move in a very similar path for the 
first second of the trajectory. However, separation has 
occurred in less than 1.3 inches of vertical displacement and 
the deviation of the motion was determined by the forces and 
torques prior to separation. 

Specific Objectives - Having examined the significant parameters 
which affect the determination of the spring/damper constants, 
the specific objectives of the study can now be identified. 

A matrix of spring and damper constants will be tested to 


12 




FIGURE e 





locate those combinations which produce non-dynamic pitching 
torque profiles. From this elimination procedure the criteria 
of trajectory deviation will be used to reduce the selection 
to a recomended set of spring/damper constants. The effects of 
dynamic pressure and frict'on coefficient will be determined 
and the results compared with similar trajectories determined 
with non-spring-damper-mass techniques and also with no active 
constraint model . 

4.0 RESULTS 

A selection of spring/damper constants was initiated by using 

6 5 

the previously described data for K = 10 and C = 10 as a start- 
ing point. A matrix of other test conditions was established to 
locate the largest values that could be used in the simulation. 

As can be seen in Figure 7 some of the cases that failed were due 
to excessively large damper constants although acceptable spring 
constants were available. Theoretically, the failure of the 

Q 

K = 10 case was predicted by the data in Figure 1. The matrix 
was completed by producing nine cases which represented a 3x3 
box of useable cases. 

Pitching Torque Profiles - The pitching torque profiles from 
these nine cases are presented in Figures 8 and 9 for the left 
and right attach, respectively. The data are grouped by damper 
constants for each of the three spring constants. Both figures 
show similarities in the frequency of the oscillation but the 
right attach, wliich remained in contact longer, shows not only 
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nxjre oscillations but larger amplitudes. As can be seen in the 
figures, each of the subplots has a spring/damper constant 
combination which produces at least one non-dynamic torque 

5 

profile. In the case of C = 10 ,on“.y the pitching torques with 

7 5 

K = 10 shows any indication of Dynamics but the K = 10 case 

suggests a sluggishness to separate. Several of the curves 

"bottom out" indicating that the ball separated and rejoined 

5 

the socket much like a bounce. The C = 10 data, unlike the 
others, shows a large initial spike in pitching torque which is 
quickly damped out. 

Three cf the nine cases are considered to have satisfied the 
criteria for non-dynamic pitching torque profiles while two 
others can be considered as close. A summary of these results 
is presented in Figure 10 where regions of degrees of accept- 
ability and failure are constructed based on the test cases. 

This figure shows an acceptable region surrounded by a probably 
acceptable region. The selection of a suitable spring/damper 
constant combination should, therefore, be made from these two 
regions. Failure cases resulted from too large of an integration 
interval for the chosen constants. 

All of the nine cases were made with typical conditions for 

dynamic pressure and friction coefficient of the ball and 

socket. The effect of varying these parameters was examined 

6 6 

for the same constants of K = 10 and C = 10 and the matrix 
of these cases are presented in Figure 11. A range of dynamic 
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pressures from 2 psf to 14 psf and friction coefficients from 
.1 to .3 are included in this matrix. The pitching torque 
profiles for the range of dynamic pressures at y = .2 is 
presented in Figures 12 and 13 and although the absolute level 
of the data has changed (since portions of the forces are 
derived from air loadu) the profile shape remains basically 
unchanged. A similar conclusion can be reached for variation 
In the friction coefficient which shows very little change. 

The largest difference occurs in the initial spike in the 
pitching torques as shown in Figures 14 and 15. 

Trajectory Profiles - Although the pitching torque profiles 
can be used to narrow down the selection, the final choice 
must consider the integration of these profiles to produce 
acceptable displacement and rotation trajectories. The dis- 
placement trajectories for the center of the ball for the left 
and right attach are presented in Figues 16 and 17. These 
trajectories cover the first inch of travel which is just prior 
to final separation. These figures show that the trajectories 
for most of the constant combinations follow a path similar to 
the geometric track of axial displacement. Two obvious exceptions 
are for K = 10^ at C = 10^ and 10^ which are two of the 
three combinations which satisfies the pitching torque criteria. 
Only one of the three cases (K = 10^ and C = 10^) shows a close 
match in the displacements when compared with the expected 
geometric track which is snown on the right edge of the Figures 
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EFFECT OF FRICTION COEFFICIENT ON PITCHING TORQUES 
OF, ORB ITER /ET AR AHACH 
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16 and 17. The results from comparing the rotation trajectories 
Is shown in Figure 18 tor both the left and right attach. These 
curves show the relative insensitivity of this trajectory to the 
extremes of the nine spring/damper constant configurations during, 
the first full second of the separation. This difference appears 
as a A0 error of about .015 degrees at the end of one second. 

Constraint Model Comparisons - An evaluation of the ABIND con- 
straint model with the other alternatives has been made in an 
earlier analysis (Reference 1). Additional comparisons are made 
here which emphasize the similarities and differences previously 
identified. These comparisons ai^e made with the constrained and 
unconstrained simulations from the MDAC-STL program as well as 
the unconstrained simulation in SVDS and are presented in Figures 
19 and 20. Figure 19, which presents the displacement trajectories 
of the ball center in the pitch plane, shows that the two con- 
straint models follow a very similar path except for the time 
difference in reaching a given point. It also demonstrates one 
of the primary reasons for using a constraint model which is to 
simulate the restricted motion of the ball in the aft direction. 
This restriction shows as a reduction of about half of the 
axial displacement by the time the ball has left the socket 
(AZ = 6 inches) . 

Figure 20 presents a somewhat different conclusion as it shows that 
the constraint models produce more rotation and less vertical dis- 
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placement than the unconstrained cases. This would be ex- 
pected as the friction coefficient converts the potential energy 
of the vertical displacement into rotating the vehicles with 
respect to each other. The constraint model also changes the ro- 
tation center from the c.g. of each vehicle to the ball/socket 
attach point. However, more importantly, it shows that use of 
a constraint simulation produces nearly 4 times as much additional 
rotation with the MDAC-STL model (A0 = .093 degrees or 38%)than 
with the ABIND model (A0 = .024 degrees or 10%). This would 
indicate that one major difference between the constraint models 
is in the ability to convert potential energy to rotation. Were 
it not for the fact that the displacement trajectories produce 
the desired reduction in axial displacement there would be little 
advantage in using the ABIf^D constraint model in the SVDS program. 




5.0 CONCLUSIONS 

An examination of the factors which modify the simulation of 
binding in the aft attach points of the Orbiter and External 
Tank during separation has been completed. The factors con- 
sidered were both internal (spring and damper constants) and 
external (friction coefficient and dynamic pressure). The 
criteria for the examination was based on an acceptable, inte- 
gration interval, a non-dynamic pitching torque-time history 
profile and a smooth displacement of the trajectories of the 
center of the ball as it rotates and translates out of the 
socket. As a result of these considerations the following 
conclusions can be stated: 

1) An 'estimate of the maximum integration interval allowed 

in the study can be made using spring-mass dynamic analyses 

2) An integration interval of .01 seconds is acceptable for 
all springs with constants less than 2.08 x 10^ Ibs/ft for 
all tank weights exceeding 80,000 lbs. 

3) Separation times for the left and right side are dependent 
on initial out-of-plane conditions. The right side is in 
contact longer for an Orbiter with a positive roll rate 
with respect to the Tank. 

4) The torque profiles are similar for either aft attach point 

5) The pitching torque is the most significant parameter for 
assessing spring/damper effectiveness. 

6) Based on analysis of the pitching torque time histories, 
only three of the test cases produce non-dynamic profiles. 


33 


7) Pitching torques show little sensitivity to variations of 
friction coefficients between .1 and .3. 

8) Pitching torque profiles have different magnitudes but 
similar shape for dynamic pressures between 2 and 14 psf. 

9) Only one of the three profile-acceptable test cases exhibited 
smooth displacements in the axial direction during separation. 

10) All nine test cases produce only slight deviations in the 
rotation trajectories. 

11) Only one of the nine test cases satisfies all of the 
established criteria (K = 10^, C = 10^) . 

12) Differences in displacement trajectories due to constraint 
model methods are small. 

13) ABIND does not produce as severe a change in the rotational 
trajectories as the MDAC-STL constraint model. 

14) A choice of spring/damper constants exist which allow ABIND 
to model the forces and torques at the aft attach points. 

6.0 RECOMM EN DATIONS 

The analyses performed herein has demonstrated that differences 
between the MDAC-STL and SVDS constraint models exist. It is 
therefore recommended that a series of cases comparing the two 
different constraint simulation models with the unconstrained 
simulation be made using initial conditions near the fonner 
left hand boundary of the separation window. Comparison of 
these cases will indicate the sensitivity of the boundary to the 
constraint model . 
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